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Due to their superfluid properties some compact astrophysical objects, like neutron or quark 
stars, may contain a significant part of their matter in the form of a Bose-Einstein Condensate. 
Observationally distinguishing between neutron/quark stars and Bose-Einstein Condensate stars 
is a major challenge for this latter theoretical model. An observational possibility of indirectly 
distinguishing Bose-Einstein Condensate stars from neutron/quark stars is through the study of 
the thin accretion disks around compact general relativistic objects. In the present paper, we 
perform a detailed comparative study of the electromagnetic and thermodynamic properties of the 
thin accretion disks around rapidly rotating Bose-Einstein Condensate stars, neutron stars and 
quark stars, respectively. Due to the differences in the exterior geometry, the thermodynamic and 
electromagnetic properties of the disks (energy flux, temperature distribution, equilibrium radiation 
spectrum and efficiency of energy conversion) are different for these classes of compact objects. 
Hence in this preliminary study we have pointed out some astrophysical signatures that may allow 
to observationally discriminate between Bose-Einstein Condensate stars and neutron/quark stars, 
respectively. 
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I. INTRODUCTION 

Since its proposal by Bose flj, and generalization by 
Einstein @, the quantum statistics of integer spin parti¬ 
cles (bosons) did represent a fundamental field of study 
in both theoretical and experimental physics. One of 
the most important property of bosonic systems is their 
phase transition to a condensed state, in which all parti¬ 
cles are in the same quantum ground state. This quan¬ 
tum bosonic system is called a Bose-Einstein Condensate 
(BEC), and from a physical point of view it is charac¬ 
terized by a sharp peak over a broader distribution in 
both coordinate and momentum space. The quantum 
explanation for this behavior is that in a BEC parti¬ 
cles become correlated with each other, and their wave¬ 
lengths overlap. Note that the correlation means that the 
thermal wavelength At is greater than the mean inter¬ 
particle distance l. It occurs at a critical temperature 
T c < 2nh 2 n 2 / 3 /mks, where m is the mass of an individ¬ 
ual condensate particle, n is the number density, and kn 
is Boltzmann’s constant M- A coherent quantum state 
in the bosonic system develops either when the particle 
density p is high enough, or when the temperature T is 
sufficiently low. From the experimental point of view, 
the Bose-Einstein Condensation can be detected by the 
observation of a sharp peak in the velocity distribution, 
which always appears when the system is below the crit¬ 
ical temperature, T < T c . 
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In the laboratory the Bose-Einstein Condensation was 
observed first in 1995 in dilute alkali gases, such as vapors 
of rubidium and sodium. To obtain the condensation the 
gases were confined in a magnetic trap, and cooled down 
to very low temperatures [9J. The experiment producing 
a laboratory BEC did represent a major achievement in 
experimental condensed matter physics, and the confir¬ 
mation of the old and important predictions in theoretical 
statistical physics of Bose and Einstein, [gf. Note that 
in recent years quantum degenerate gases have been cre¬ 
ated by a multitude of experimental methods, including 
combination of laser and evaporative cooling techniques. 
Thus the observation of the Bose-Einstein Condensation 
did open several new lines of multidisciplinary research 
at the border of atomic, statistical and condensed matter 
physics M- 

The Bose-Einstein condensation processes are assumed 
to play an important role in the understanding of many 
fundamental processes in condensed matter physics. For 
example, superfluidity of low temperature liquids, like 
3 He, can be explained by assuming a Bose-Einstein Con¬ 
densation process It is interesting to note that ex¬ 
perimental observations as well as quantum theoretical 
calculations estimate the condensate fraction no at T = 0 
for superfluid helium to be only around no ~ 0.10, a very 
small amount. Hence, since a strongly correlated pair of 
fermions behaves approximately like a boson, the aris¬ 
ing liquid helium superfluidity can be interpreted as re¬ 
sulting from the Bose-Einstein Condensation of coupled 
fermions. A very similar physical model can be used to 
describe the transition to a superconducting state in a 
solid material (metal). Superconductivity can thus be 
interpreted as the condensation of electrons (or holes) 
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into Cooper pairs. The formation of the Cooper pairs 
drastically reduces the resistance caused by the motion 
of electrons in metals, and thus leads to the formation of 
a superconductor [gj. 

Since Bose-Einstein condensation is a phenomenon 
that has been observed and thoroughly studied in terres¬ 
trial laboratories, the possibility that it may also occur in 
bosonic systems existing on astrophysical or cosmological 
scales cannot be rejected a priori. Thus, it was proposed 
that dark matter, which is required to explain the dynam¬ 
ics of the neutral hydrogen clouds at large distances from 
the centers of the galaxies, and which is assumed to be a 
cold bosonic gravitationally bounded system, could also 
exist in the form of a Bose-Einstein Condensate 0. a 
systematic study of the properties of Bose-Einstein Con¬ 
densed galactic dark matter halos was initiated in a 
and the astrophysical and cosmological implications of 
the existence of Bose-Einstein Condensed dark matter 
have been investigated in detail recently 0 - 

From a theoretical point of view it was shown in El 
that by introducing the Madelung representation of the 
wave function, the dynamics of the Bose Einstein Con¬ 
densate dark matter halo is described by the continu¬ 
ity equation, and the hydrodynamic Euler equations of 
the standard classical fluid mechanics. Hence, Bose- 
Einstein Condensed dark matter can be described the¬ 
oretically as a gas, with the pressure and density re¬ 
lated by a barotropic equation of state. In the case of 
self-interacting condensate dark matter with quartic self¬ 
interaction potential, the equation of state is polytropic 
with index n = 1 11 ill. 


Bose-Einstein Condensation could play an important 
role in nuclear and quark matter physics, in the frame¬ 
work of the so-called Bardeen-Cooper-Schrieffer (BCS) 
to Bose-Einstein Condensate (BEC) crossover. From a 
theoretical point of view it is expected that at ultra-high 
densities nuclear matter exists in the form of a degener¬ 
ate Fermi gas of quarks. The Cooper pairs of quarks 
form near the Fermi surface a Bose-Einstein Conden¬ 
sate. Hence high density nuclear matter represents, from 
a physical point of view, a so called color superconductor 
13]. When the attractive interaction between fermions is 
strong enough, and the temperature drops below the crit¬ 
ical temperature, the fermions condense into the bosonic 
zero mode, and form a quark BEC 14]. As a first step to¬ 
wards the formation of the BEC the fermions must form 
a BCS state, which can be realized when the attractive 
interaction between particles is weak. This system ex¬ 
hibits superfluid properties, which are characterized by 
the existence of an energy gap for single particle excita¬ 
tions. The energy gap is created by the formation of the 
Cooper pairs. On the other hand a BEC is formed when 
the attractive interaction between fermions is extremely 
strong. This interaction first leads to the formation of 
bound particles (bosons), which at some critical temper¬ 
ature T c start to condense into the bosonic zero mode. It 
is important to mention that the BCS and BEC states are 
smoothly connected (crossover), without a phase transi¬ 


tion between the two phases |15|. For a recent review of 
the BCS-BEC crossover see fl6j . 


The possibility of the existence of some forms of Bose- 
Einstein Condensates in neutron stars has been consid¬ 
ered a long time ago (see Glendenning 0 for a detailed 
discussion). One possibility for the formation of a BEC in 
a dense neutron star is the condensation of the negatively 
charged mesons, leading to the replacement of electrons 
with very high Fermi momenta by mesons 0. Bose- 
Einstein Condensation of kaons and anti-kaons in com¬ 
pact objects was also investigated in detail 0 , 0 . It 
turns out that the presence of pion or kaon condensates 
may have at least two important effects on the global 
properties of dense neutron stars. Firstly, Bose-Einstein 
Condensation softens the equation of state (EOS) of the 
stellar matter above the critical density for the onset of 
condensation. An important consequence of the soften¬ 
ing of the EOS is the reduction of the maximum neutron 
star mass. On the other hand, due to the softening of 
the EOS, the central density increases significantly. Sec¬ 
ondly, the condensation of mesons would lead to con¬ 
siderably enhanced neutrino luminosities, much higher 
than those of normal neutron matter. The increase in 
neutrino luminosity has important consequences on the 
neutron star cooling 0. Another particle, which may be 
present inside neutron stars, and which may form a Bose- 
Einstein Condensate, is the H-dibaryon. The H-dibaryon 
is a doubly strange six-quark composite. It has zero spin 
and isospin, and a baryon number B = 2 0. Neutron 
star matter may also contain an important fraction of 
A hyperons, neutral subatomic hadrons. They consist of 
one up, one down and one strange quark, and they are la¬ 
belled A 0 0| . The A hyperons may also combine to form 
H-dibaryons (2lj| . Thus, H-matter Bose-Einstein conden¬ 
sates may also exist at the center of very dense neutron 
stars, where the matter density is extremely high 0 - 
An interesting possibility, is that neutrino superfluidity, 
as suggested by Kapusta 0, may also lead to Bose- 
Einstein Condensation inside neutron stars 0 ] . 


Thus theoretical results in nuclear matter physics in¬ 
dicate that the possibility of the existence of some forms 
of Bose-Einstein Condensed matter inside compact astro- 
physical objects, like neutron or quark stars, or even the 
existence of stars formed entirely from a pure BEC, can¬ 
not be excluded a priori. The properties of pure BEC 
stars have been considered in [24]. It was shown that 
stars formed of Bose-Einstein Condensates with parti¬ 
cle masses of the order of two neutron masses (Cooper 
pair) and scattering length of the order of 10-20 frn have 
maximum masses of the order of 2 M 0 , maximum cen¬ 
tral densities of the order of 0.1 — 0.3 x 10 16 g/cm 3 and 
minimum radii in the range of 10-20 km. Hence Bose- 
Einstein Condensed stars can form a large class of stable 
astrophysical objects, whose basic astrophysical param¬ 
eters (mass and radius) sensitively depend on the mass 
of the condensed particle, and on the scattering length. 
It was also suggested that the recently observed neutron 
stars with masses in the range of 2 — 2.4Mg (Vela X- 
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1, 4U 1700-377, and the black widow pulsar B1957+20) 
are Bose-Einstein Condensate stars. Further properties 
of BEC stars have been considered in [§5|. 

The structure of static and rotating BEC stars was 
investigated numerically in ,261] by solving the Gross- 
Pitaevskii-Poisson system of coupled differential equa¬ 
tions. It was shown, with longer simulation runs, that 
within the computational limits of the simulation the 
BEC stars are stable. The physical properties of the self- 
gravitating Bose-Einstein Condensate were investigated 
in both non-rotating and rotating cases. 

A fundamental theoretical problem is to find some 
clear astrophysical signatures, going beyond the global 
stellar parameters (mass and radius) that could indicate 
the presence of a BEC inside a neutron star, or to give a 
firm observational evidence for the existence of pure BEC 
stars. It is the purpose of the present paper to suggest 
that such a specific signature, indicating the presence of 
a pure BEC star, may indeed exist, and it can be ob¬ 
tained from the study of the radiation emission of thin 
accretion disks that usually form around compact general 
relativistic objects. 

The growth of most astrophysical objects is determined 
by mass accretion. Due to the presence of interstellar 
matter, accretion disks are generally formed around com¬ 
pact objects. Accretion disks are well known observa- 
tionally, representing flattened astronomical structures. 
They are made of rapidly rotating hot gas, slowly spiral¬ 
ing onto a central massive and dense object. The grav¬ 
itational energy of the gas motion is a source of heat, 
generated by the internal stresses and the dynamical fric¬ 
tion. A small fraction of the heat is converted into ra¬ 
diation, which partially escapes. The radiation emission 
cools down the accretion disk. Therefore important in¬ 
formation about the accretion disk physics comes from 
the radiation emitted from the disk. The radiation, de¬ 
tected in the radio, optical or X-ray frequency bands, 
allows astronomers to analyze its electromagnetic spec¬ 
trum, and its time variability. Thus essential results 
about the physics of the disks can be obtained from ob¬ 
servations. 

The cooling of the disk via the electromagnetic ra¬ 
diation emission from its surface represents an efficient 
mechanism that prevents the extreme heating of the disk. 
On the other hand, the thermodynamic equilibrium es¬ 
tablished in this way allows the disk to stabilize its thin 
vertical size. Usually the inner edge of the thin disk is 
located at the marginally stable orbit of the compact ob¬ 
ject potential. Hence in higher orbits the hot gas has a 
Keplerian motion USUI. Since the electromagnetic ra¬ 
diation emission from the thin disk is determined by the 
external gravitational potentials, which in turn are de¬ 
termined by the EOS of the dense neutron or quark mat¬ 
ter in the star, astrophysical observations of the emission 
spectra from accretion disks may lead to the possibility of 
directly testing the equation of state of the dense matter 
inside compact general relativistic objects. 

The emissivity properties of the accretion disks have 


been used to investigated, and obtain distinctive as¬ 
trophysical signatures, for large classes of compact as¬ 
trophysical objects, including naked singularities 29], 
gravastars [3C| and wormholes a Specific electro¬ 
magnetic disk signatures in different modified gravity 
theories, such as f(R) gravity, brane world models, 
Chern-Simons models and the Horava-Lifshitz theory 
were considered j32| , while the properties of accretion 
disks around rotating and non-rotating neutron, quark, 
boson or fermion stars have been analyzed in [33H35| . 

In the present paper we perform a comparative study 
of the rotational properties and of the disk emission prop¬ 
erties for five neutron star EOSs, and of the quark matter 
bag model EOS with the polytropic n = 1 BEC equation 
of state. The main goal of this preliminary investigation 
is to point out towards the possible existence of some ob¬ 
servational signatures that may distinguish between these 
different classes of compact objects. We begin our analy¬ 
sis by considering the global rotational properties of the 
considered compact objects. To obtain the equilibrium 
configurations of the rotating neutron, quark and BEC 
stars we use the Rotating Neutron Star (RNS) code, as 
introduced in [36], and discussed in detail in [37]. The 
RNS software provides the metric potentials for different 
types of equations of state of compact rotating general 
relativistic objects. We investigate three different cases, 
corresponding to stellar models with fixed mass and an¬ 
gular velocity, stellar models rotating at Keplerian fre¬ 
quencies, and stellar models with fixed central density 
and fixed polar radius to equatorial radius ratio. As a 
next step in our study we use the exterior metrics to ob¬ 
tain the physical properties of the accretion disks for the 
considered equations of state. Particular signatures ap¬ 
pear in the electromagnetic spectrum of the BEC stars, 
thus leading to the theoretical possibility of directly test¬ 
ing, and discriminating, the BEC equation of state of the 
dense matter by using astrophysical observations of the 
emission spectra from accretion disks. 

The present paper is organized as follows. In SectionllD 
we briefly review the basics of the Bose-Einstein Con¬ 
densation and discuss the properties of Newtonian BEC 
stars. The electromagnetic and thermodynamic proper¬ 
ties of accretion disks around compact general relativis¬ 
tic objects are described in Section ITTT1 The equations of 
state of the neutron, quark and BEC matter, as well as 
the global astrophysical properties of the considered stel¬ 
lar models are obtained and discussed in Section [TV] The 
electromagnetic and thermodynamic properties (flux, lu¬ 
minosity and temperature distribution) of the accretion 
disks around neutron, quark and BEC stars are obtained 
in Section El We discuss our results and conclude our 
study in Section ED 


II. BOSE EINSTEIN CONDENSATION 

The most important characteristic of a quantum sys¬ 
tem of N interacting bosons in the condensed state is 
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that most of the particles lie in the same single-particle 
quantum state. For a bosonic quantum system consist¬ 
ing of an extremely large number of particles, the cal¬ 
culation of the physical parameters of the ground state 
with the direct use of the Hamiltonian is generally very 
difficult, due to the high computational cost. A signifi¬ 
cant simplification of the many-body type computations 
can be achieved with the use of some approximate meth¬ 
ods. One such approximate semi-analytic method is the 
mean field analysis of the quantum condensate. The ba¬ 
sic idea of the mean field method is the separation of 
the Bose-Einstein Condensate contribution to the total 
bosonic field operator. In the following in our analysis we 
assume that the bosonic astrophysical system consists of 
scalar (zero spin) particles with non-zero mass. From a 
simplifying physical point of view we assume that when 
the system makes a transition to a Bose-Einstein Con¬ 
densed phase, the range of Van der Waals-type scalar 
mediated interactions among particles becomes infinite. 


A. The Gross-Pitaevskii equation 


In the second quantization approach the Hamiltonian 
describing a many body system of interacting bosons, 
confined by an external potential V ex t, is given by 



where $ (fi) and <f> + (fi) are the bosonic field operators, 
and V (fi — fi) is the two-body interatomic potential, re¬ 
spectively i-1- The bosonic field operators annihilate 
and create a particle at the position fi. V ro t (fi) is the 
potential associated to the rotation of the condensate. 

In order to obtain a significant simplification of the 
mathematical formalism a number of approximate meth¬ 
ods have been developed. One such approximate and 
simplifying semi-analytic approach is the mean field de¬ 
scription of the condensate. In this approach the con¬ 
densate contribution to the bosonic field operator is sep¬ 
arated out. For a uniform gas confined in a volume V , 
BEC occurs in the single particle state <I>o = 1\/V, having 
zero momentum. The field operator can then be decom¬ 
posed in the form <h (fi) = \JN/V + <3E>' (fi). By treating 
the operator (fi) as a small perturbation, the first order 
theory for the excitations of the interacting Bose gases 
can be fully developed I'l I* . 

In the general case of a non-uniform and time- 
dependent configuration the Heisenberg representation of 
the field operator is given by 

$ (r, t) = if (fi, t) + (fi, t) , (2) 


where if(r,t), also called the condensate wave function, 
is the expectation value of the field operator, if (fi, t) = 


^ (fi, t) j. It is important to note that the wave function 

if (fi, t) is a classical field, and its absolute value fixes 
the particle number density of the condensate through 
the relation p(fi,t) = \if (fit )| 2 . On the other hand the 
normalization condition for the condensate wave function 
is N = J p(fi t) d 3 fi where N is the total number of 
particles in the Bose-Einstein Condensate. 

The condensate wave function satisfies the following 
equation of motion, which can be obtained from Heisen¬ 
berg equation corresponding to the many-body Hamilto¬ 
nian given by Eq. ©, 


ih— $ (r, t) = &,H 

dt K ' 1 L ’ . 


— X-V 2 + V ro t (fi) + V ex t (fi) + 

2m 


dr 4> + (fi, t)V (fi — fi) $ (fi, t) 


$ (fi, t). 


( 3 ) 


The zeroth-order approximation to the Heisenberg 
equation is obtained by replacing $ (r, t) with the con¬ 
densate wave function if. For short distances this is in 
general a poor approximation for computing the integral 
containing the particle-particle interaction V (fi — fi). 
However, it is important to note that in a dilute and 
cold gas, at low energy only binary collisions are im¬ 
portant. These collisions can be characterized, indepen¬ 
dently of the exact form of the two-body potential, by a 
single physical parameter, the s-wave scattering length a. 
Therefore, one can obtain a very good approximation by 
replacing the generally unknown potential V (fi — fi) with 
an effective interaction potential V (fi — fi) = XS (fi — fi), 
with the coupling constant A determined by the scatter¬ 
ing length l a via the relation A = 47r h 2 a/m, where m is 
the mass of the condensate particles. With the use of the 
effective potential the integral in the bracket of Eq. © 
can be easily calculated to give A \if (fi, t)\ 2 . Hence the 
resulting approximate equation of motion for the conden¬ 
sate is the Schrodinger equation with a quartic nonlin¬ 
ear term 0-1 , also called the Gross-Pitaevskii equation. 
However, in order to obtain a more general and realistic 
description of the Bose-Einstein Condensate stars, one 
may also assume an arbitrary non-linear self-interaction 
term of the form g (\if (fi,t) | 2 j = g (p) i- 

Therefore the non-relativistic generalized Gross- 
Pitaevskii equation, describing a gravitationally trapped 
rotating Bose-Einstein condensate, is given by 


ih—if(r,t) = 


- ^V 2 + V rot (fi + V^t (fi) + 

2m 


i (IV’(g^)I 2 ) 


if ( fi ., t), 


( 4 ) 


where we have denoted g' = dg/dp. As for V ext ( fi ), we 
assume that it is the gravitational potential, denoted for 
simplicity by V, V ex t = V. The gravitational potential 
satisfies the Poisson equation of Newtonian gravity 

V 2 P = 4tt Gp m , 


( 5 ) 
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where p m = mp = m \ip (r, t)\ 2 is the mass density inside 
the Bose-Einstein condensate. 


B. The hydrodynamical representation 

The physical as well as the astropliysical properties of 
a Bose-Einstein Condensate, described by the general¬ 
ized Gross-Pitaevskii equation Eq. 0 can be analyzed, 
computed and understood much more easily by using the 
so-called Madelung representation of the wave function 
P]. In this representation we write the condensate wave 
function t/j in the specific form 


ip (r, t) = \Jp (r, t) exp 




( 6 ) 


where the function S (r, t ) has the physical dimensions 
of an action. By substituting the above expression of 
(r, t) into Eq. 0, it follows that the generalized Gross- 
Pitaevskii equation decouples into a system of two first 
order partial differential equations for the two real func¬ 
tions p m and v, respectively, given by 


& Pm 

~df 


+ V ■ ( p m V ) = 0, 


(7) 


Pm 


dv 

m 


+ (v ■ V) v 



Pm^ 




W Q , 


( 8 ) 


where we have introduced the quantum potential Vq , de¬ 
fined as, 


exact form of the equation of state of the BEC depends 
on the non-linearity term g. 

In the important case of a gravitationally bounded 
Bose-Einstein Condensate with a very large number of 
particles the quantum pressure is important only near 
the external boundary of the quantum system. Hence 
in all astrophysically important situations the quantum 
pressure term is much smaller than the non-linear inter¬ 
action term Ell¬ 
in the most studied approach to the Bose-Einstein 
Condensates, the non-linearity term g is quadratic, and 
is given by 

9(Pm) = ^ M 4 = yPm> ( 12 ) 

where uq = Airh 2 l a /m i-0- The corresponding equa¬ 
tion of state of the condensate is 

P(p) = Kp 2 . (13) 

Therefore we have obtained the very important result 
that the equation of state of the standard Bose-Einstein 
Condensate with quartic non-linearity is a polytrope, 
with index n = 1. 

In the following for simplicity we consider only the 
case of the Bose - Einstein condensates with quartic non¬ 
linearity. In this particular case the physical properties 
of the condensates are also relatively well known from nu¬ 
merous laboratory experiments. From theoretical point 
of view it is important to note that their properties can 
be described in terms of only two free parameters, the 
mass to of the condensate particle, and the scattering 
length a, describing the particle interaction, respectively. 


C. Masses and radii of Newtonian static BEC stars 


fl 2 V 2 
2 TO 


(9) 


the velocity of the Bose-Einstein Condensed fluid v, given 

by 


TO 


( 10 ) 


and we have denoted 

p(Ell) = g>(eil)eH-g(Ell). ( 11 ) 

V TO / V TO / TO V TO / 

From its definition we can immediately see that the 
velocity field v is irrotational. Thus it satisfies the condi¬ 
tion V x v = 0. Therefore we obtain the important result 
that in the Madelung representation the equations of mo¬ 
tion of the Bose-Einstein Condensate in a gravitational 
field take the form of the equation of continuity, and of 
the hydrodynamic Euler equations, respectively. Hence 
a Bose-Einstein Condensate in an external gravitational 
field can be described as a gas whose density and pressure 
are related by a barotropic equation of state Q-Q. The 


Since the EOS of a cold BEC star is a polytrope of 
index n = 1, all the physical properties of the star can 
be derived from the well-known Lane-Emden equation, 
given by [24j 


1 d ( d9 \ 

v W 


= - 0 , 


(14) 


where 8 is a dimensionless variable defined via p = p c 9 , 
with p c the central density of the star, and £ is a dimen¬ 
sionless radial coordinate. The general non-singular solu¬ 
tion of Eq. (m is well known, and can be represented as 
9 (£) = sin£/£. We define the dimensionless radius of the 
star by the condition 9 (£i) = 0, giving = n. Hence the 
physical radius R of the Bose-Einstein Condensate star 
can be represented as a function of the catering length 
and the fundamental physical constants as FFlj 


R = TT 




1/2 


1 fm / 


2 m r 


-3/2 


km. (15) 


For to = 2 m n and a = 1 fm the numerical value of the 
radius of the BEC star is R ~ 7 km. It is interesting 
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to note that the radius of the gravitationally bounded 
pure Bose-Einstein Condensate star is independent on 
the central density and on the mass of the star, and is 
fully determined by the physical characteristics a and m 
of the condensate only. 

The mass of the BEC star can be found as [Hj 

M =n<^) »• < ie > 


stable circular orbits in stationary and axially symmetric 
geometries. The general formalism for obtaining the con¬ 
stants of the motion (energy and angular momentum), 
the effective gravitational potential and the radii of the 
marginally stable orbits is presented in some detail. Then 
we briefly review the basic equations describing the elec¬ 
tromagnetic radiation emission from accretion disks, and 
present the general expressions for the flux, temperature 
distribution and luminosity of the disk. 


or, equivalently, 


M = 1.84 


10 16 g cm -3 


a \ 3 / 2 / m \ 9//2 
1 fm/ \2m n ) 


Mq. 


(17) 

For m = 2m n , a = 1 fm and p c = 5 x 10 15 g/cm 3 , the 
mass of the condensate is M ft 0.92 Mq. On the other 
hand by taking for the mass of the condensate particle 
the numerical value m = m* K = m n / 10, where m* K is the 
kaon mass, we obtain an extremely high maximum mass 
of M = 63.50M 0 for the kaon condensate star, with a 
corresponding radius of R = 355 km. Note that the mass 
of the static condensate can be expressed in terms of the 
radius and central density by 


M = - p c R 3 . 

7r 


(18) 


The above important equation shows that the mean den¬ 
sity of the star p = 3M/4irR 3 can be obtained from the 
central density of the condensate by the simple relation 
P = 3p c /7T 2 . 

Under a scaling of the parameters m, a and p c so that 
m — > aim , a —> Pc — > a^p Cl the radii and the masses 

ofthe Newtonian BEC stars have the scaling properties 


R -> oq 3/2 a 1 / 2 R, M -A ctq 9/2 « 2 /2 a 3 M. (19) 

It is important to point out that general relativistic 
effects impose strong constraints on the structure and 
global parameters of the BEC stars, and that the values 
obtained by the simple Newtonian estimates may highly 
exceed the stability limit imposed by the general rela¬ 
tivistic analysis. 


III. ELECTROMAGNETIC RADIATION 
PROPERTIES OF THIN ACCRETION DISKS IN 
STATIONARY AXISYMMETRIC SPACETIMES 


In order to investigate the accretion disk properties 
around neutron, quark and BEC stars we first briefly in¬ 
troduce the general formalism that allows us to describe 
the electromagnetic radiation properties of thin accre¬ 
tion disks in stationary axisymmetric spacetimes. In our 
presentation we closely follow the approach developed in 
[33! . As a first step in our study we consider the ba¬ 
sic kinematic properties of massive particles moving in 


A. Massive particle motion in stationary and 
axially symmetric spacetimes 


In this work our main emphasis is on the analysis of the 
physical properties and characteristics of particles that 
form a thin accretion disk, and move in circular stable 
orbits around general relativistic compact objects. The 
exterior geometry created by a central dense object is 
assumed to be stationary and axially symmetric, given 
in full generality by the following metric, 

ds 2 = g tt dt 2 + 2 g t(/) dtdcj) + g rr dr 2 + gee dd 2 + d(f) 2 . 

( 20 ) 

In the equatorial approximation, i.e., \9 — 7r| <C 1, which 
we adopt in our study, all the metric functions gtt, gt</>, 
Qrr, gee and g^ depend on the radial coordinate r only 
[27, 28]. In the following we denote the square root of 
the determinant of the metric tensor by \J—g. 

In order to determine the electromagnetic properties of 
the disk we first obtain the radial dependence of the angu¬ 
lar velocity II, of the specific energy E and of the specific 
angular momentum L for particles moving in circular or¬ 
bits around compact objects in the geometry given by 
Eq. (l20l) . All these physical parameters can be obtained 
from the geodesic equations, which for the considered 
metric take the following form [29|-(32j 


dt 

ds 

dcj) 

ds 



Eg^rf, + Lgtcf, 

9tp-9ttgM ’ 

Egtt/> + Lgu 

9tct> ~ 9ttgp<j> 

1 A' 2 ff00 + 2 ELg tt f, + L 2 g tt 

9%-9tt9^ 


( 21 ) 

( 22 ) 

(23) 


We define an effective potential term V e ff(r) as 


Veff{r) =-l + 


E 9$$ dr 2ELgt ( j > -{- L gtt 

9t(p ~ 9ttg4>4> 


(24) 


For particles moving in stable circular orbits in the 
equatorial plane the potential V e ff(r) must satisfy the 
following two important conditions: V e ff(r) = 0 and 
V e ff, r (r) = 0, respectively, where a coma denotes the 
derivative with respect to the radial coordinate r. From 
these conditions we obtain the specific energy, the specific 
angular momentum and the angular velocity of particles 
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moving in circular orbits for the case of spinning general 
relativistic compact objects in the form [33| 


E 

L 

n 


gtt + gt<pQ 

\J gtt g(/n/>Q^ 

gt<p t 

\J gtt 2gt(pQ g^^Q? 

d(f) _ 9t(j),r d" \/( 9t(f>,r9tt,r9cl>(l>,r 

dt 9cf>cl),r 


(25) 

(26) 
(27) 


The marginally stable orbit around the central object 
can be determined from the supplementary condition 
V e ff. rr(f) = 0, which provides the following important 
relationship 


E g(/)(j),rr T 2E Lgt^^rr -\- L gtt,rr ( 9t(ft 9tt9(j>(l>) 7 rr — 0- 


(28) 


Once the metric coefficients gtt, gt 4 > and g ^ are ex¬ 
plicitly given, by inserting Eqs. (E5l) - (E71) into Eq. (051) . 
and numerically solving this equation for the radial co¬ 
ordinate r, we obtain the radii of the marginally stable 
orbits of massive particles in stable circular orbits around 
a rotating general relativistic high density compact ob¬ 
ject p9|-[33|. 


B. Electromagnetic emissivity of thin accretion 
disks 


four-velocity an have a rotational velocity Q = d(j)/dt , 
a specific energy E, and a specific angular momentum L. 
In the steady state thin disk model all these quantities 
depend only on the radii r of the orbits. The particles 
form a disk of an averaged surface density E, which is 
obtained as the average of the rest mass density po of 
the gas in the vertical direction. The matter in the disk 
is described by an anisotropic fluid source. The den¬ 
sity po of the disk, the energy flow vector q M and the 
dissipative part of the energy-momentum tensor are 
all measured in the averaged rest-frame. In this specific 
frame the specific heat and heat transfer processes are 
neglected. Then, one important parameter characteriz¬ 
ing the disk structure is the surface density of the disk 


27, 


E(r) = [ {po)dz, (29) 

J-H 

which is obtained as the averaged rest mass density (po) 
over At and 27 t. Another important disk parameter is 
the torque 



representing the average of the component (fj) of the 
energy-momentum tensor over At and 27 t. The time and 
orbital average of the energy flux vector gives the radia¬ 
tion flux F(r) from the disk surface as 


In the following we concentrate on thin accretion disks, 
that is, accretion disks having their vertical size negligible 
as compared to their horizontal extensions. This means 
that the disk height H, which we define as the maxi¬ 
mum half thickness of the disk in the vertical direction, 
is much smaller than the characteristic radius R of the 
disk, H <C R- We further assume that the thin disk is 
in hydrodynamical equilibrium. Moreover, we fully ne¬ 
glect the possible effects of the pressure gradient, and 
of the vertical entropy gradient in the disk. The heat 
generated by internal stresses and dynamical friction is 
emitted over the disk surface, and this emission prevents 
the disk from heating up at extremely high temperatures. 
On the other hand, this equilibrium determines the disk 
to stabilize its thin vertical size. An important astro- 
physical parameter, the inner edge of the thin disk, is 
located at the marginally stable orbit of the central ob¬ 
ject gravitational potential. Hence the accreting matter 
has a Keplerian motion in higher orbits. 

In the steady state accretion disk model, the mass ac¬ 
cretion rate Mq is a constant that does not change in 
time. An effective physical description of the disk proper¬ 
ties is obtained by averaging all the quantities describing 
the orbiting gas over a characteristic time scale, e.g. At, 
over the azimuthal ang le A</> = 2-7T, and over the height 
H, respectively fill. |38L [39]. 

As we have already seen, the particles moving in Ke¬ 
plerian orbit around the dense compact object with the 


Fir) = (9 Z )- (31) 

The energy-momentum tensor of the matter in the disk 
is represented in its standard form according to 

= poit^u" + 2 + t^, (32) 

where the four-velocity of the matter satisfies the con¬ 
ditions = 0, u^ 1 ' = 0. The four-vectors of 

the energy and angular momentum flux are defined by 
-E* = Ttf(d/dt) v and J * = TJf (d / d4>Y , respectively. 
Note that the structure equations of the thin disk can be 
derived by integrating the conservation laws of the rest 
mass, of the energy, and of the angular momentum of 
the gas forming the disk, respectively [13, [391 • From the 
equation of the rest mass conservation, V ^(pqu^) = 0, it 
follows first the important result that the time averaged 
mass accretion rate is independent of the disk radius, 

Mq = —2tt ^/~^gT,u r = constant. (33) 

The conservation law = 0 of the energy can be 

reformulated in the integral form, more interesting from 
a physical point of view, as [331 

[M 0 E - 27ry^f2W/] , r = 4 t t^^FE. (34) 

The above equation shows that the energy transported 
by the mass flow, MqE, and the energy transported by 









the dynamical stresses in the disk, 2tt • v / = gf 2 W ( £ r , is ex¬ 
actly compensated by the energy radiated away from the 
surface of the disk, 4 n^/—gFE. The law of the angu¬ 
lar momentum conservation, = 0, indicates the 

equilibrium between the three forms of the angular mo¬ 
mentum transport, 

[M 0 L - 2tt rW/], r = An^gFL . (35) 

By eliminating r from Eqs. (El and (1551) . with 
the use of the universal energy-angular momentum re¬ 
lation dE = QdJ for circular geodesic orbits, written 
in the form E r = £lL^ ri the flux F of the electromag¬ 
netic e nerg y emitted by the disk surface is obtained as 

Fir) = - , _ n,r ~ [ (E-QL)L r dr , (36) 

(E - UL ) 2 J rin 

where ?y n is the inner edge of the disk. In our study we 
assume that ri n = r ms . Note that the flux depends of the 
specific energy, angular momentum and angular velocity 
of the gas motion around the central general relativistic 
compact object. 

The gas forming the accretion disk in the steady-state 
thin disk model is assumed to be in thermodynamical 
equilibrium. Therefore the radiation emitted by the disk 
surface may be considered as a perfect black body radia¬ 
tion, depending on the temperature only, with the energy 
flux given by 

F(r) = a SB T\r ), (37) 


where asB is the Stefan-Boltzmann constant. Once the 
radiative flux is known, from Eq. (1371) we obtain the tem¬ 
perature distribution on the disk surface. 

The observed luminosity L (y) of the disk surface has 
a redshifted black body spectrum, given by |34j, 


L (y) = And 2 ! (b) = 


8nh cos 7 


r r f 


B^rd(f>dr 


I o exp ( hue/ksT ) - 1 ’ 


(38) 


In Eq. (1551) h is Planck’s constant, k B is Boltzmmann’s 
constant, d is the distance to the source (disk), 1(b) is the 
Planck distribution function, B e is the frequency of the 
emitted radiation, 7 is the disk inclination angle, and 
and rf indicate the positions of the inner and outer edge 
of the disk, respectively. In the natural system of units 
with h = c = k B = 1 , we obtain for the disk luminosity 
the expression 


r T f r 

L (b) = 16n 2 cos 7 / / 

Jr in Jo 


B^rd(j)dr 


exp (27 TB e /T) — 1 ’ 


(39) 


In the following we take the upper limit of integra¬ 
tion in Eq. (1381) as infinity, that is, we assume 77 —> 00 . 
Moreover, we expect that the flux from the disk surface 


vanishes at r —> 00 . This condition is independent on 
the geometry of the general relativistic compact object. 
The frequency of the emitted radiation is given by the 
relation B e = b( 1 + z), where the redshift factor can be 
written as 

1 + Sir sin d> sin 7 

1 + z=—= v -L - . (40) 

v 9tt U 2 ^^ 

In Eq. (l40l) we have neglected the light bending [4(1 . 

This approximation works well for small inclination an¬ 
gles, but it is not as good for large inclination angles 
(edge-on disks). 

Another important characteristics of the mass accre¬ 
tion processes is the efficiency e with which the central 
object converts the mass of the gas into radiation. The 
efficiency e is defined as the ratio of the rate of the energy 
of photons emitted from the disk surface, and the rate at 
which mass-energy is transported to the central high den¬ 
sity general relativistic object [13, [Hj. Both energies are 
measured at infinity. If all photons escape to infinity, the 
efficiency depends on the specific energy measured at the 
marginally stable orbit r ms only, 

e = 1 — E ms . (41) 

For Schwarzschild black holes the efficiency e is about 
6 %. This result is independent on whether the photon 
capture by the black hole is considered, or not. For 
rapidly rotating black holes for which the capture of ra¬ 
diation by the black hole is ignored, eis equal to 42%. 
For a Kerr black hole with the photon capture explicitly 
considered the efficiency is 40% [28| . 

Note that the fluxes and the emission spectra of the ac¬ 
cretion disks around compact objects satisfy some simple, 
but important, scaling relations. Such scaling relations 
can be found when one considers the scaling transforma¬ 
tion of the radial coordinate, given by r —> r = r/M , 
where M is the mass of the central object. Under such 
a scaling transformation in general the metric tensor co¬ 
efficients are invariant, while the specific energy, the an¬ 
gular momentum and the angular velocity transform as 
E —> E, L —>• ML and fl —> tt/M, respectively. The flux 
scales as F(r) —> F(r)/M 4 , giving the simple scaling law 
of the temperature as T(r) T (?) /M. Another im¬ 
portant rescaling is the transformation of the frequency 
of the emitted radiation as v —> b = b/M. Under this 
rescaling of the radiation frequency the luminosity of the 
disk scales as L(b) —>■ L (b) /M |33j]. On the other hand, 
since the flux is proportional to the accretion rate Mo, 
we obtain the important result that an increase in the 
accretion rate leads to a linear increase of the radiation 
flux from the disk. 

From an observational point of view, the use of the 
thermal radiation from thin disks to test the nature of 
a compact general relativistic object is often called the 
continuum-fitting method, and was proposed in [43[ to 
investigate the observational consequences of black hole 
spin in X-ray binaries. The standard thin accretion disk 
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model introduced in j27t was used to analyze the emis- 
sivity of accretion disks around black holes. In the case 
of Kerr black holes the observed disk spectrum, how¬ 
ever, bears several important corrections to the simple 
formula given by Eq. (l36l) . Since the X-rays are emitted 
on the hot inner disk, in this region electron scattering 
may dominate over the free-free absorption. As a con¬ 
sequence the color temperature may be greater than the 
effective temperature, with the inner disk radiating like 
a diluted black body. The general relativistic corrections 
near the event horizon (gravitational redshift and focus¬ 
ing) can cause the observed color temperature and the 
integrated flux to be different from the true local values. 
The continuum-fitting method was extensively used to 
observationally study the spin of black holes j44l |. 


IV. EQUATIONS OF STATE AND PHYSICAL 
PARAMETERS OF THE NEUTRON, QUARK 
AND BEC STARS 

To set the stage for our study, in the present Section 
we introduce first the equations of state of the dense nu¬ 
clear matter considered in our comparative study of the 
neutron, quark and Bose-Einstein Condensate stars, and 
present some of the basic dynamical properties of the cor¬ 
responding rotating stellar models. In order to compare 
the properties of the BEC stars with other compact ob¬ 
jects we choose five equations of state of neutron matter, 
and the bag model equation of state for quark matter. 
We restrict our study to rotating BEC stars satisfying 
the n = 1 polytropic equation of state. The astrophysical 
parameters of the stars described by the adopted equa¬ 
tions of state are also presented in detail. In order to 
compare the properties of the stellar models we consider 
three classes of models. In the first model we fix the val¬ 
ues of the mass and of the angular velocity for all stars. 
As a second case we consider compact stars rotating at 
the maximum (Keplerian) frequencies. And finally, as a 
third class of models we consider stars having the same 
central density and fixed ratio of the equatorial and polar 
radii, respectively. 


A. Equations of state of the neutron, quark and 
BEC matter 

In order to obtain a consistent and realistic physical de¬ 
scription of the rapidly rotating general relativistic com¬ 
pact stars, as a first step we have to adopt the equations 
of state for the dense neutron, quark, and Bose-Einstein 
Condensate stellar matter, respectively. In the present 
comparative study of the physical properties of the ac¬ 
cretion disks onto rapidly rotating compact general rela¬ 
tivistic objects we adopt the following equations of state 

a „ 

1) Akmal-Pandharipande-Ravenhall (APR) EOS |45j ]. 
As a first example of an equation of state of dense 


neutron matter we consider EOS APR, which was ob¬ 
tained numerically by using the variational chain sum¬ 
mation methods and the Argonne Vis two-nucleon in¬ 
teraction, respectively. Note that boost corrections to 
the two-nucleon interaction, giving the leading relativis¬ 
tic effect of order (v/c) 2 , as well as three-nucleon inter¬ 
actions, are also included in the basic nuclear Hamil¬ 
tonian. The dense matter density range described by 
EOS APR is from 2 x 10 14 g/cm 3 to 2.6 x 10 15 g/cm 3 . 
The maximum mass limit in the static case for this 
EOS is 2.20 Mq. To obtain a full description of the 
neutron star properties the APR EOS is joined to the 
composite Baym-Bethe-Pethick (BBP) (e/c 2 > 4.3 x 
10 11 g/cm 3 ) [i|| - Baym-Pethick-Sutherland (BPS) (10 4 
g/cm 3 < 4.3 x 10 n g/cm 3 ) [47| - Feynman-Metropolis- 
Teller (FMT) (e/c 2 < 10 4 g/cm 3 ) [H} equations of state, 
respectively. 

2) Douchin-Haensel (DH) EOS [49}. EOS DH repre¬ 
sents a complete equation of state of the neutron star 
matter. It describes both the neutron star crust as well 
as its liquid core. It is constructed by using the effective 
nuclear interaction SLy of the Skyrme type. The cor¬ 
responding interaction potential describes very well the 
properties of very rich neutron matter. On the other 
hand the structure of the crust, and its EOS, is obtained 
in the zero temperature approximation only, and under 
the assumption of the ground state composition. The 
EOS of the liquid core is calculated assuming (minimal) 
tipe/u composition. The density range is from 3.49 x 10 11 
g/cm 3 to 4.04 x 10 15 g/cm 3 . The minimum and maxi¬ 
mum masses of the static neutron stars for the DH EOS 
are O.O94M 0 and 2.05 Mq, respectively. 

3) Shen-Toki-Oyamatsu-Sumiyoshi (STOS) EOS [si)} . 
The STOS equation of state of nuclear matter is obtained 
by using the relativistic mean field theory with nonlinear 
a and to terms. It can be used in a wide neutron matter 
density and temperature range, with various proton frac¬ 
tions. The EOS was specifically designed for the use of 
supernova explosion simulation and for the calculations 
of the neutron star properties. To compute the proper¬ 
ties of inhomogeneous nuclear matter, where heavy nuclei 
are formed together with free nucleon gas, the Thomas- 
Fermi approximation is used. In the present paper we 
consider only the zero temperature STOS EOS, namely 
the EOS for T = 0, and we denote it as STOSO. 

4) Relativistic Mean Field (RMF) equations of state 
with isovector scalar mean field in the presence of the 
5-meson held- RMF soft and RMF stiff EOSs [5l}. Even 
that in symmetric nuclear matter the 5-meson mean held 
vanishes, this held still can significantly inhuence the 
properties of asymmetric nuclear matter in neutron stars. 
Note that the Relativistic Mean Field contribution due 
to the 5-held to the nuclear symmetry energy is smaller 
than zero. Therefore the energy per particle of neutron 
matter is larger at high densities, as compared to the 
case when no (5-held included. Also, the proton fraction 
of /3-stable matter increases. Splitting of proton and neu¬ 
tron effective masses due to the 5-held can significantly 
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affect transport properties of neutron star matter. These 
equations of state can be parameterized by the coupling 
parameters C 2 = gl/m 2 a , Cl = g 2 u /m 2 u , b = b/gl and 
c = c/gl, respectively, where m„ and m u are the masses 
of the respective mesons, and b and c are the coefficients 
in the potential energy U (a) of the cr-field. The soft RMF 
EOS is parameterized by Cl = 1.582 fm 2 , C 2 = 1.019 
fm 2 , b = —0.7188 and c = 6.563, while the stiff RMF 
EOS is parameterized by C 2 = 11.25 fm 2 , C 2 = 6.483 
fm 2 , b = 0.003825 and c = 3.5 x 10 -6 , respectively. 

5) Baldo-Bombaci-Burgio (BBB) EOS [52|. The BBB 
EOS is an EOS for asymmetric nuclear matter. It is 
derived from the Brueckner-Bethe-Goldstone many-body 
theory with explicit three-body forces taken into ac¬ 
count. Two EOS’s are obtained, one corresponding to 
the Argonne AY14 (BBBAV14), and the other to the 
Paris two-body nuclear force (BBBParis), implemented 
by the Urbana model for the three-body force. The max¬ 
imum static mass configurations are M max = 1.8 Mq and 
M max = 1.94Mq when the AV14 and Paris interactions 
are used, respectively. The onset of direct Urea processes 
occurs at particle number densities n > 0.65 fm -3 for 
the AV14 potential, and n > 0.54 fm -3 for the Paris po¬ 
tential. The comparison with other microscopic models 
for the EOS shows noticeable differences, which can also 
influence significantly the neutron star properties. The 
density range for this EOS is from 1.35 x 10 14 g/cm 3 to 
3.507 x 10 15 g/cm 3 . 

6 ) Bag model equation of state for quark matter - Q 
EOS [53l |. For the description of the quark matter we 
adopt a simple phenomenological description, based on 
the MIT bag model equation of state. Hence we assume 
that in quark matter the pressure p is related to the en¬ 
ergy density p by 


P=^(p-4B)c 2 , (42) 

where the parameter B , called the bag constant, is 
the difference between the energy density of the pertur¬ 
bative and non-perturbative Quantum Chromodynamic 
vacuum. For the bag constant we adopt the numerical 
value B = 4.2 x 10 14 g/cm 3 [53j . 

7) The Bose-Einstein Condensate equation of state 
(BEC) EOS. For the Bose-Einstein Condensate stars we 
adopt the n = 1 polytropic equation of state, given by 

P(p) = Kp 2 , (43) 


with 


K = 


2n h 2 a 
m 3 


= 0.1856 x 10 5 



(44) 


where m n = 1.6749 x 10 ~ 24 g is the mass of the neu¬ 
tron. The behavior of the BEC EOS, as well as of the 
corresponding stellar models essentially depends on the 
ratio a/m 3 of the scattering length and of the conden¬ 
sate particle mass. In the present study we restrict our 
analysis to three values of the ratio (a /1 fm) ( m/2m n ) 3 : 


(a/1 fm) (m/2m n ) 3 = 15, (a/1 fm) (m/2m n ) 3 = 30, 
and (a/1 fm) (m/2m ri )~ 3 = 50, respectively. We denote 
the corresponding equations of state of the BEC matter 
by BEC15, BEC30, and BEC50, respectively. 

The variation of the considered equations of state with 
respect to the density of the dense stellar matter is rep¬ 
resented in Fig. |I] 



P [g cnU 


FIG. 1: Variation of the pressure as a function of the density 
for the considered equations of state of the neutron, quark 
and BEC dense matter. 

The low density behavior of the BEC15 and BEC30 
EOSs is significantly different as compared to the other 
considered equations of state, indicating the possibility 
of the existence of stable BEC stars at pressures lower 
than the neutron matter pressure. In the adopted range 
of densities the BEC EOSs show a linear pressure-density 
dependence. 


B. Astrophysical parameters of the neutron, quark 
and BEC stars 

In quasi-isotropic coordinates the metric outside a ro¬ 
tating compact general relativistic star can be repre¬ 
sented as [36( 

ds 2 = -e^ +p dt 2 + e 2a (df 2 + r 2 d0 2 ) + e^~ p f 2 sin 2 9 x 
(d(j> — udt) 2, , (45) 

where the metric potentials 7 , /?, a and the angular ve¬ 
locity of the stellar fluid ui, measured relative to the lo¬ 
cal inertial frame, are all functions of the quasi-isotropic 
radial coordinate f, and of the polar angle 6. The 
RNS code computes numerically the metric functions 
in a quasi-spheroidal coordinate system, as functions of 
the parameter s = r/(r + f e ), where f e is the equa¬ 
torial radius of the star, which we have converted into 
Schwarzschild-type coordinates r according to the equa¬ 
tion r = r exp [(7 — p) /2]. To obtain the radius of the 
marginally (or innermost) stable circular orbits r ms we 
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use a truncated form of the analytical approximation 
given as 


MI, 


^ « l- 0.54433 g - 0.22619 9 * 2 + 0.17989Q 2 - 
6M 

0.23002q 3 + 0.26296 9 Q 2 - 0.29693q 4 + 
0.44546<? 2 Q 2 , (46) 

where q = J/M 2 and Qi = —Af 2 /Af 3 , respectively, and 
where J is the spin angular momentum, and Af 2 is the 
quadrupole moment. 


mass stable BEC star is of the same order as the radii of 
the ’’standard” stars, R e « 8 — 15 km. Significant differ¬ 
ences appear in the maximum masses of the stars, with 
the maximum allowable mass of the BEC50 star varying 
between 2.8Af 0 (static case) and 3.0M 0 , for the rapidly 
rotating star with r p /r e = 0.7. For the BEC15 star, 
the maximum mass of the stable configuration ranges be¬ 
tween 1.5 -Mq (static case), and 1.6M 0 , for r p /r e = 0.7. 
Hence a specific mass-radius relation provides a first dis¬ 
tinctive feature of the Bose-Einstein Condensate stars, as 
compared to the considered neutron and quark stars. 


1. Mass-radius relation for neutron, quark and BEC stars 

The mass-radius relation of the compact general ob¬ 
jects for the considered equations of state of the dense 
nuclear matter, M = M ( R e ), where R e is the circumfer¬ 
ential radius at the equator, are presented, for four fixed 
values of the ratio r p /r e , where r p is the polar radius, 
and r e is the equatorial radius of the star, in Figs. [2] and 
[3l respectively. 

As one can see from the Figures, the BEC stars form 
a distinct class of stellar objects, as compared to the 
group of ’’standard” neutron and quark stars. The mass- 
radius relation is systematically shifted to the right re¬ 
gion of the M = M (R e ) relation, has a specific shape, 
and indicates a much larger radius for the BEC star. 
The radius increases with increasing a/m 3 , so that for 
(a/1 fm) (m/2m n )~ 3 = 50 and r p /r e = 0.7 the radius of 
the maximum mass stable BEC star is of the order of 30 
km. Thus a first distinctive signature of rapidly rotat¬ 
ing BEC stars is their bigger radius, as compared to the 
’’standard” neutron and quark stars, indicating a large 
value of the coefficient a/m 3 . On the other hand, for 
(a/1 fm) ( m/2m n ) < 20, the radius of the maximum 

I 


The BEC15 EOS does not allow stellar masses of the or¬ 
der of 1.8 Mq. The mass corresponding to an angular 
velocity of Q w 5 x 10 3 s _1 is 1.541 Mq, with a small 
equatorial radius of the order of 6 km, and a high central 
density p c = 5 x 10 15 g/cm 3 . Configurations with 1.8M 0 
can be obtained for the BEC30 and BEC50 EOSs. They 
have equatorial radii of the order of 22 and 28 km. A sim¬ 
ilar, but still smaller equatorial radius r e = 18 km, can 
be found only for the STOSO EOS. The BEC30, BEC50 
and STOSO configurations have similar central densities. 
The BEC30 and BEC50 stars have the highest moment 


2. Models with fixed mass and angular velocity 

In analyzing the emissivity properties of the accretion 
disks we consider three types of stellar models, whose 
main astrophysical properties are presented in the fol¬ 
lowing in a tabular form. In all Tables, p c is the central 
density, M is the gravitational mass, Mq is the rest mass, 
R e is the circumferential radius at the equator, H is the 
angular velocity, Q p is the angular velocity of a particle 
in circular orbit at the equator, T/W is the rotational- 
gravitational energy ratio, cJ/GMq is the angular mo¬ 
mentum, / is the moment of inertia, 4> 2 gives the mass 
quadrupole moment M 2 so that Af 2 = c 4 3> 2 /(G 2 Af 3 AfJ;), 
h + is the height from the surface of the last stable co¬ 
rotating circular orbit in the equatorial plane, h- is the 
height from surface of the last stable counter-rotating cir¬ 
cular orbit in the equatorial plane, w c /H is the ratio of 
the central value of the potential ui to f2, r e is the coordi¬ 
nate equatorial radius, and r p /r e is the axes ratio (polar 
to equatorial), respectively. 

The physical properties of the neutron, quark and BEC 
stars with fixed mass, M ss 1.8Af 0 and angular velocity 
H«5x 10 3 s _1 are presented in Table [1] 


of inertia /, and highest angular momentum per unit 

cJ/GM , with EOS STOSO having the closest values of 
these parameters. 

3. Models rotating at Keplerian frequencies 


As a second astrophysical model we consider the case of 
the neutron, quark and BEC stars rotating at Keplerian 
frequencies. The physical properties of this class of stars 
are presented in Table m 


The most massive stable star for this model, with mass 
AT « 3.5 M e is obtained for the STOSO EOS. The BEC50 


star has a comparable, but somewhat lower mass, of 
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FIG. 2: Mass-circumferential radius relation for the rotating neutron, quark and BEC stars for r p /r e = 1 (left figure) and 
r p /r e — 0.9 (right figure). 




FIG. 3: Mass-circumferential radius relation for the rotating neutron, quark and BEC stars for r p /r e = 0.8 (left figure) and 
r p /r e = 0.7 (right figure). 


EOS 

APR 

BBBAV14 

BBBParis 

BEC15 

BEC30 

BEC50 

DH 

Q 

RMFsoft 

RMFstiff 

STOSO 

p c [10 15 g/cm 3 ] 

1.20 

2.15 

1.70 

5.00 

0.50 

0.19 

1.29 

0.895 

2.00 

0.57 

0.369 

M[Mq\ 

1.798 

1.805 

1.799 

1.541 

1.811 

1.8 

1.805 

1.806 

1.525 

1.803 

1.854 

Mo [Mo] 

2.083 

2.108 

2.095 

1.735 

1.944 

1.895 

2.054 

1.953 

1.705 

2.003 

2.016 

R e [km] 

11.437 

10.35 

10.706 

8.35 

24.917 

30.661 

12.02 

11.979 

11.344 

15.812 

21.038 

SI[10 3 s -1 ] 

5.008 

5.034 

5.039 

5.04 

4.269 

2.581 

5.008 

5.006 

4.999 

5.009 

4.902 

SI P [10 3 s” 1 ] 

12.473 

14.414 

13.705 

18.341 

3.983 

2.961 

11.593 

11.781 

11.674 

7.972 

5.374 

T/W[10 -2 ] 

3.449 

2.375 

2.712 

1.113 

9.756 

9.429 

3.46 

4.836 

3.312 

8.981 

14.911 

cJ/GMq 

1.142 

0.953 

1.01 

0.487 

2.233 

2.418 

1.155 

1.389 

0.826 

2.011 

2.953 

J[10 45 g cm 2 ] 

2.005 

1.663 

1.761 

0.849 

4.596 

8.233 

2.028 

2.438 

1.452 

3.528 

5.294 

<P 2 [10 43 g cm 2 ] 

8.263 

4.325 

5.365 

0.978 

69.829 

143.288 

8.433 

14.988 

6.717 

43.932 

106.248 

h+ [km] 

0.00 

3.39 

2.93 

3.82 

0.00 

0.00 

1.65 

0.00 

0.00 

0.00 

0.00 

h- [km] 

7.958 

8.352 

8.166 

6.845 

1.296 

0.00 

7.47 

0.00 

0.00 

0.00 

7.859 

cu c /S2[10 _1 ] 

5.809 

6.667 

6.343 

7.769 

4.159 

3.025 

5.856 

5.252 

5.519 

4.526 

4.103 

r e [km] 

8.521 

7.408 

7.781 

5.841 

22.087 

27.848 

9.11 

9.024 

8.909 

12.872 

18.013 

r p /r e 

0.894 

0.925 

0.915 

0.958 

0.54 

0.623 

0.885 

0.86 

0.883 

0.728 

0.54 


TABLE I: Astrophysical parameters of the neutron, quark and BEC stars with mass M « 1.8Mq, rotating at an angular 
velocity S2 « 5 x 10 3 s^ 1 . 
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EOS 

APR 

BBBAV14 

BBBParis 

BEC15 

BEC30 

BEC50 

DH 

Q 

RMFsoft 

RMFstiff 

STOSO 

p c [10 15 g/cm 3 ] 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

M[Mq] 

1.93 

1.646 

1.671 

1.286 

2.311 

3.18 

1.784 

2.807 

1.841 

2.786 

3.492 

M 0 [M q \ 

2.183 

1.83 

1.86 

1.381 

2.557 

3.567 

1.982 

3.316 

2.039 

3.229 

4.147 

Re [km] 

15.493 

15.279 

15.408 

16.831 

20.132 

22.291 

16.249 

17.152 

17.793 

18.435 

19.148 

fl[10 3 s- 1 ] 

8.398 

7.948 

7.919 

6.07 

6.105 

6.059 

7.439 

9.031 

6.912 

7.746 

7.991 

fl p [10 3 s” 1 ] 

8.405 

7.948 

7.919 

6.07 

6.105 

6.06 

7.439 

9.031 

6.912 

7.746 

7.991 

T/W[ 10” 2 ] 

15.064 

14.074 

14.347 

9.932 

9.71 

9.458 

11.49 

21.989 

16.232 

15.493 

14.885 

cJ/GMq 

2.88 

2.081 

2.168 

1.138 

3.295 

6.038 

2.154 

7.183 

2.913 

5.948 

8.975 

1[10 45 g cm 2 ] 

3.014 

2.301 

2.406 

1.648 

4.744 

8.757 

2.545 

6.99 

3.704 

6.748 

9.871 

< I? 2 [ 10 43 g cm 2 ] 

47.829 

37.726 

40.021 

25.417 

54.725 

84.065 

33.907 

130.893 

73.347 

99.348 

123.857 

h+ [km] 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.87 

1.63 

h- [km] 

10.41 

7.397 

7.695 

1.893 

8.805 

15.843 

6.918 

21.205 

9.532 

17.732 

24.495 

Wc/flflO" 1 ] 

5.789 

5.176 

5.203 

4.171 

6.047 

7.332 

5.439 

6.978 

5.154 

6.905 

7.945 

r e [km] 

12.29 

12.597 

12.679 

14.81 

16.422 

17.093 

13.379 

12.044 

14.728 

13.672 

13.056 

r p /r e 

0.52 

0.527 

0.524 

0.568 

0.574 

0.578 

0.558 

0.447 

0.494 

0.522 

0.543 


TABLE II: Astrophysical properties of the neutron, quark and BEC stars rotating at Keplerian frequencies, for p c = 10 15 
g/cm 3 . 


around 3.2M 0 . The RMFstifF and the quark equation 
of state Q have Keplerian masses of the order of 2.8 M@, 
greater than the Keplerian mass of the BEC15 EOS, 
M = 2.31M©. The smallest Keplerian mass is obtained 
for the BEC15 EOS. The BEC stars have the biggest 
equatorial radii, ranging between 15 and 17 km. For the 
neutron and quark star EOSs the highest equatorial ra¬ 
dius is obtained for the STOSO EOS, but which is still 
smaller than the equatorial radius of the BEC15 model. 
The equatorial radius of the RMSsoft EOS almost coin¬ 
cides with the Keplerian equatorial radius of the BEC15 
EOS, with the RMSsoft EOS having a much greater mass 
than the BEC15 EOS one. The Keplerian model of the 


Similarly to the previous cases, the higher mass for 
stars having fixed central density and a polar to equato¬ 
rial radius ratio r p /r e = 0.85 is obtained for the STOSO 
EOS, M = 2.9792A7©. The mass of the BEC50 EOS 
model is lower, with a numerical value of around 2.9M©. 
The lowest mass value is obtained for the BEC15 EOS, 
M = 1.11M©. The mass of the BEC30 EOS model, 
M = 2.068M© is significantly smaller than the masses 
obtained for the Q EOS, with M = 1.924M©, and 
for the RMFstiff EOS, having a corresponding mass of 
M = 2.266 Mq. The BEC stars have the larger equato¬ 
rial radii, of the order of 13-18 km, while the DH EOS 
has an equatorial radius r e ~ 12.6 km, the largest equa¬ 
torial radius in this class of stellar models. The BEC15 
EOS has the smallest moment of inertia and angular mo¬ 
mentum per unit mass, with the highest values of these 
parameters obtained for STOSO EOS. However, the an¬ 


STOSO EOS has the greatest moment of inertia, and an¬ 
gular momentum per unit mass, with the BEC15 having 
the smallest values for these quantities. 


4- Models with fixed central density and r p /r e = 0.85 

The astrophysical parameters of the neutron, quark 
and BEC stars with a fixed central density p c = 10 15 
g/cm 3 , and a fixed ratio of the polar to the equatorial 
radius r p /r e are presented in Table m 


gular momentum per unit mass of the BEC50 EOS just 
slightly exceeds the angular momentum per unit mass 
of the RMSstiff EOS, but it is smaller than the cJ/GM 
value of the quark EOS Q. 


V. ELECTROMAGNETIC AND 
THERMODYNAMIC SIGNATURES OF 
ACCRETION DISKS AROUND BOSE-EINSTEIN 
CONDENSATE STARS 

In the present Section we consider the electromagnetic 
signatures of the accretion disks around neutron, quark 
and BEC stars. We consider a comparative study in¬ 
volving three distinct classes of stellar models. The first 
model corresponds to accretion disks formed around com¬ 
pact general relativistic objects with fixed masses, of the 
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EOS 

APR 

BBBAV14 

BBBParis 

BEC15 

BEC30 

BEC50 

DH 

Q 

RMFsoft 

RMFstiff 

STOSO 

Pc[10 15 g/cm 3 ] 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

M[Mq\ 

1.539 

1.322 

1.334 

1.116 

2.068 

2.907 

1.536 

1.924 

1.393 

2.266 

2.979 

Mo [Mq\ 

1.735 

1.464 

1.478 

1.196 

2.286 

3.264 

1.702 

2.265 

1.537 

2.626 

3.557 

Re [km] 

11.723 

11.507 

11.569 

12.825 

15.803 

17.899 

12.619 

12.102 

12.686 

14.145 

15.395 

f![10 3 s -1 ] 

5.273 

5.030 

4.995 

4.132 

4.201 

4.221 

4.911 

5.359 

4.385 

4.951 

5.207 

O p [10 3 s" 1 ] 

11.247 

10.761 

10.725 

8.417 

8.263 

8.034 

10.036 

11.974 

9.576 

10.195 

10.153 

T/W[10~ 2 ] 

4.702 

4.517 

4.536 

3.808 

3.917 

3.996 

4.286 

5.245 

4.51 

4.879 

5.211 

o 

C$ 

s; 

©to 

1.006 

0.747 

0.762 

0.517 

1.643 

3.239 

0.96 

1.628 

0.844 

2.169 

3.86 

7[10 45 g cm 2 ] 

1.677 

1.305 

1.341 

1.100 

3.436 

6.745 

1.718 

2.670 

1.693 

3.850 

6.518 

4>2[10 43 g cm 2 ] 

10.785 

8.946 

9.236 

7.664 

17.988 

29.527 

10.39 

16.897 

12.332 

21.675 

32.598 

h+ [km] 

0.00 

0.00 

0.00 

0.00 

0.00 

3.286 

0.00 

0.00 

0.00 

0.00 

5.403 

h- [km] 

5.896 

3.947 

4.034 

0.689 

7.105 

13.618 

4.835 

0.00 

0.00 

11.145 

17.492 

Wc/fffKT 1 ] 

4.954 

4.419 

4.427 

3.755 

5.634 

6.998 

4.875 

5.613 

4.3 

6.135 

7.366 

r e [km] 

9.25 

9.405 

9.447 

11.08 

12.507 

13.173 

10.177 

8.93 

10.476 

10.433 

10.387 

r p /r e 

0.85 

0.85 

0.85 

0.85 

0.85 

0.85 

0.85 

0.85 

0.85 

0.85 

0.85 


TABLE III: Astrophysical parameters for neutron, quark and BEC stars for p c = 10 15 g/cm 3 and r p /r e = 0.85. 


order of M « 1.8M 0 , rotating at an angular speed of 
Q « 5 x 10 3 s 1 . The second case corresponds to accretion 
disks formed around stars rotating at the maximal Keple- 
rian frequency. And, finally, we also consider a third class 
of models, in which the accretion disk is located around 
a star with fixed central density and polar to equatorial 
radius ratio. For all these three cases we consider the 
disk emissivity properties, which are strongly dependent 
of the equation of state of the dense matter inside the 
star. 


A. Electromagnetic spectrum from accretion disks 
around rotating neutron, quark and BEC stars with 
fixed mass and angular velocity 

We begin our analysis of the electromagnetic signa¬ 
tures of accretion disks around compact general rela¬ 
tivistic objects by considering the case of stars with 
fixed mass, M « 1.8M 0 , and angular velocity of 11 ~ 
5 x 10 3 s _1 . The variations of the electromagnetic flux, 
disk temperature and luminosity for neutron, quark and 
BEC stars with fixed mass and angular velocity are pre¬ 
sented in Figs. Ml 

As one can see from the Figures, the BEC stars form 
a distinct group with respect to the neutron and quark 
stars included in the study. The flux emitted by the ac¬ 
cretion disks, presented in Fig. [H is the smallest for the 
BEC50 and BEC30 equations of state, respectively. The 
inner disk edge for the BEC50 and BEC30 EOSs is lo¬ 
cated at r/M ~ 12 and r/M ~ 9, respectively, at the 
highest distance from the central object for all consid¬ 
ered stars. The maximum value of the flux is obtained 
for EOS RMFsoft, and this maximum flux value is about 
twelve and five times bigger than the maximum flux val¬ 
ues from the BEC50 and BEC30 EOSs. The flux emitted 



FIG. 4: Electromagnetic fluxes from accretion disks gravitat¬ 
ing around compact general relativistic objects, having differ¬ 
ent equations of state, with mass M = 1.8 M®, and rotating 
at an angular velocity of fl = 5 x 10 3 s -1 . 


by the star with EOS BEC15 has high values of the flux, 
comparable with those from EOS RMFsoft, but this can 
be explained by the lower mass of the star, 1.541M 0 , and 
the corresponding scaling of the flux. However, for the 
BEC15 EOS the inner edge of the disk is located at a 
distance of around 5.5 x r/M from the central object. 

The temperature distribution in the disk, shown in 
Fig. M generally follows the same distribution as for the 
flux profiles, with the BEC stars having some specific 
distinctive features. The lowest maximum disk tempera¬ 
ture is obtained for the BEC50 and BEC30 EOSs, with 
the maximum temperature located at around 17 x r/M 
and 13 x r/M , respectively. For most of the neutron and 
quark stars the maximum disk temperature is reached at 
r/M ft: 8. Some specific distinctive features also appear 
for EOS STOSO, with a maximum temperature located 

















15 




FIG. 5: Temperature distribution of the accretion disks 
around compact general relativistic objects, having different 
equations of state, with mass M = 1.8 Mq, and rotating at 
an angular velocity of Q. == 5 x 10 3 s . 


FIG. 7: Electromagnetic fluxes from accretion disks around 
compact general relativistic objects with different equations 
state rotating at Keplerian frequencies. 
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FIG. 8: Temperature distribution of the accretion disks 
around compact general relativistic objects with different 
equations of state rotating at Keplerian frequencies. 


FIG. 6: Luminosity of the accretion disks around compact 
general relativistic objects, having different equations of state, 
with mass M = 1.8 Mq, and rotating at an angular velocity 
of O = 5 x 10 3 s _1 . 


B. Electromagnetic signatures of accretion disks 
gravitating around neutron, quark and BEC stars at 
Keplerian frequencies 


at around 11 x r/M, relatively close to the temperature 
maximum for EOS BEC30. The highest temperature of 
the disk is reached by EOS RMFsoft. 

The maximum value of the luminosity vL{y) of the 
disk, presented in Fig. [Gl is blue shifted for the BEC50 
and BEC30 EOSs, and it is reached at a smaller fre¬ 
quency. The maximum values of the luminosity are 
smaller for these BEC EOSs by a factor of around 1.8. 
Hence a shift in the position of the luminosity maximum, 
and a lower value of the luminosity give two specific sig¬ 
natures that could help identify BEC 50 and BEC30 stars 
via the study of the luminosity of accretion disks around 
stars with known physical parameters. 


As a second example of specific electromagnetic sig¬ 
natures from accretion disks around neutron, quark and 
BEC stars we consider the case of accretion disks formed 
around compact stars rotating at Keplerian frequencies. 
The corresponding electromagnetic fluxes, the disk tem¬ 
perature distribution, and the luminosities are presented 
in Figs, midi 

From the point of view of the maximum of the flux 
emission from accretion disks around neutron, quark and 
BEC stars rotating at Keplerian frequencies, the stellar 
models considered in the present study can be roughly 
divided in two classes: disks with high flux values, and 
disks with low flux values. As can be seen from Fig. [7] 
EOSs APR, BBBAV14, Q and DH have relatively simi¬ 
lar maximum flux values, which exceeds the maximum 
flux values of EOSs RMFstiff, RMFsoft, STOSO, and 
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v[Hz] 


FIG. 9: Luminosity of the accretion disks around compact 
general relativistic objects with different equations of state 
rotating at Keplerian frequencies. 


BEC50, BEC30, and BEC15, respectively. Hence from 
the point of view of the flux maximum the BEC equa¬ 
tions of state present some similarities with the RMFs- 
tiff, RMFsoft and STOSO equations of state. However, 
important differences do appear in the localization of the 
position of the maximum of the flux. While for most of 
the equations of state the maximum is located in the re¬ 
gion 7 — 10 x r/M , the position of the maximum flux for 
EOS BEC15 is at around 13 x r/M. For this equation 
of state the inner edge of the disk is located at 9 x r/M, 
while the inner edges of the BEC50 and BEC30 EOSs 
are located in positions similar to the other neutron and 
quark matter equations of state. Hence the determina¬ 
tion of the position of the flux maximum from an accre¬ 
tion disk could give a clear indication about the nature 
of the central compact object. The temperature distri¬ 
bution of the disks around stars rotating at Keplerian 
frequencies, shown in Fig. [51 also indicates the existences 
of two distinct types of stars, with the BEC stars be¬ 
longing to the lower temperature group. However, the 
position of the temperature maximum gives a very clear 
signature of the nature of the equation of state of the 
compact object. The position of the temperature max¬ 
imum shifts towards higher distances from the central 
object with the decrease of the parameter a/m 3 . 

The luminosity of the disk, plotted in Fig. [9l shows 
that the maximum of the luminosity is reached at about 
the same frequency for all considered equations of state. 
However, the luminosity maximum has different values 
for different EOSs, with the BEC15 EOS having the 
smallest value. The BEC50 and BEC30 EOSs have simi¬ 
lar luminosity values as for the other considered neutron 
and quark matter equations of state. 



r/M 

FIG. 10: Electromagnetic fluxes from accretion disks around 
compact general relativistic objects with fixed central density 
p c = 10 15 g/cm 3 and r p /r e = 0.85. 



FIG. 11: Temperature distribution of the accretion disks 
around compact general relativistic objects with fixed central 
density p c = 10 15 g/cm 3 and r p /r e = 0.85. 


C. Electromagnetic signatures of accretion disks 
around neutron, quark and BEC stars with fixed 
central density and r p /r e = 0.85 

Finally, we consider the electromagnetic properties of 
the accretion disks around neutron, quark and BEC stars 
with fixed central density and ratio of the polar to the 
equatorial radius r p /r e = 0.85. All the stellar models 
have high angular speeds Cl. The variation with respect 
to r/M of the emitted fluxes, the temperature distribu¬ 
tion of the accretion disks, and their luminosities are pre¬ 
sented in Figs, rmirnn 

From the point of view of the electromagnetic flux dis¬ 
tribution, presented in Fig. [TUI the fluxes can be clas¬ 
sified in three groups, having high flux values (EOSs Q 
and APR), medium values (EOSs RMFstiff, RMFsoft, 
STOSO, DH, BBBAV14, BBBParis), and low flux values, 
with all the BEC EOSs belonging to this latter group. It 
is interesting to note that for all BEC EOSs the max¬ 
imum value of the flux has (approximately) the same 
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EOS 

APR 

BBBAV14 

BBBParis 

BEC15 

BEC30 

BEC50 

DH 

Q 

RMFsoft 

RMFstiff 

STOSO 

Kepler n n [km] 

e 

15.495 

0.0729 

15.280 

0.0659 

15.408 

0.0662 

16.831 

0.0498 

20.132 

0.0695 

22.291 

0.0796 

16.251 

0.0675 

17.171 

0.0807 

17.801 

0.0596 

18.438 

0.0817 

19.148 

0.0905 

= 0.85 Tin [km] 

e 

19.067 

0.0527 

16.330 

0.0524 

16.445 

0.0527 

14.030 

0.0512 

26.987 

0.0503 

37.927 

0.0505 

19.482 

0.0514 

23.327 

0.0539 

17.028 

0.0528 

28.148 

0.0527 

36.294 

0.0538 

11 = 5 x 10 3 s -1 Tin [km] 
e 

12.171 

0.0663 

13.003 

0.0670 

12.715 

0.0671 

11.952 

0.0646 

24.917 

0.0479 

30.661 

0.0389 

12.204 

0.0658 

11.979 

0.0635 

11.344 

0.0668 

15.812 

0.0635 

21.038 

0.0537 


TABLE IV: The radius of the inner disk edge r; n , and the efficiency e of the matter-radiation conversion for stars in Keplerian 
rotation, r v /r e = 0.85, and = 5 x 10 3 s^ 1 . 



v[Hz] 

FIG. 12: Luminosity of the accretion disks around compact 
general relativistic objects with fixed central density p c = 10 15 
g/cm 3 and r p /r e = 0.85. 

value. However, these maxima are located at different 
r/M, with the maximum of the BEC15 EOS positioned 
at around 15 x r/M from the central object. The maxi¬ 
mum value of the flux of the BEC stars is around three 
to four times lower than from the stars in the first group, 
with maximum fluxes (EOSs Q and APR), and around 
two times smaller as compared to the maximum flux val¬ 
ues for the EOSs of the second group. The temperature 
distribution, shown in Fig. [TT| follows the same pattern 
as the flux, with the maximum of the disk temperature 
situated at around 18 x r/M for EOS BEC15. The maxi¬ 


The first two lines of Table IIVI contain the conversion 
efficiency of the compact neutron, quark and BEC stars 
rotating at Keplerian frequencies. In this case the values 
of e for the BEC stars are in the range of 5% to 8%, 
having values comparable to the e values for neutron and 
quark stars. The highest efficiency, 9%, is obtained for 
the STOSO EOS, while the BEC50 EOS has a conversion 
efficiency close to that of the quark stars. The smallest 
e value is obtained for the BEC15 EOS, showing that 


mum temperature value is very similar for all three BEC 
equations of state. The luminosities of the disks, plot¬ 
ted in Fig. 1121 shows that the maximum of the function 
i/L(u) is reached at roughly the same frequency for all 
considered EOSs. However, there is a clear theoretical 
difference in the absolute value of the luminosity max¬ 
imum, with the BEC15 EOS having the smallest lumi¬ 
nosity. Even that EOSs BEC 50 and BEC30 have disk 
luminosities comparable with some neutron star models, 
they still belong to the class of low luminosity disks. 


D. Efficiency of radiation emission from accretion 
disks around neutron, quark and BEC stars 

An important observable physical parameter of the 
disk, which could help to observationally distinguish be¬ 
tween different classes of neutron, quark and BEC stars, 
is the efficiency e of the conversion of the accreting mass 
into radiation, given by Eq. m- The numerical values 
of e show the efficiency of the energy generating mecha¬ 
nism of the mass accretion [3^]. From a physical point 
of view, the binding energy E ms , or E e , represents the 
amount of energy released by the matter leaving the 
marginally stable orbit, or the inner edge of the disk, 
touching the surface of the star, and being transferred to 
the star. The radii of the inner disk edges for different 
equation of state and the efficiency of the radiation emis¬ 
sion is presented, for the three stellar models considered 
in the present paper, in Table Hvl 


these stars are less efficient engines for the conversion of 
the accreted mass into outgoing radiation. 

The values of e for the compact general relativistic ob¬ 
jects with fixed central density and r p /r e are given in the 
second two lines of Table [W] The values of e are slightly 
lower, as compared to the Keplerian rotation case, with 
the highest value for e obtained for the STOSO EOS. The 
smallest value of e is found for the BEC15 EOS. However, 
in this case the efficiency of the BEC50 EOS is signifi- 
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cantly smaller than the efficiency of the Q EOS, being 
slightly higher than e for APR EOS. 

The last two lines in Table IIVI show the numerical val¬ 
ues of e for disks around compact objects with mass M = 
1.8 Mq, rotating with an angular velocity of 0 ~ 5 x 10' 3 
s -1 . While the efficiency of the radiation conversion is 
almost the same for EOSs APR, BBBAV14, BBBParis, 
DH, Q, RMFsoft, and RMFstiff, it has some lower val¬ 
ues for the BEC15 EOS, and much lower values for EOSs 
BEC30 and BEC50. Relatively low values have been also 
obtained for the STOSO EOS. The efficiencies of 3.89% 
obtained for the BEC15 and BEC50 stars are the lowest 
accretion disk efficiencies obtained in the present study. 
The inner edges of the accretion disks around BEC stars 
are located at a much bigger distance from the central 
object as compared to the other classes of neutron and 
quark stars, reaching a value of 31 km for the BEC50 
EOS. Such a far away located inner edge does explain 
the low efficiency of the corresponding accretion disk. 


VI. DISCUSSIONS AND FINAL REMARKS 


The possible existence of some forms of Bose-Einstein 
Condensates in compact general relativistic objects, or 
the existence of pure Bose-Einstein Condensate stars, 
represents an intriguing, and interesting, possibility, for 
which a lot of theoretical evidence has been provided. 
From an observational point of view the most important 
differences between standard neutron or quark stars are 
represented by the differences in mass and radius. If 
the masses and radii of the compact general relativis¬ 
tic object could be measured with high accuracy, these 
measurements would put very strong direct constraints 
on the equation of state of the dense star. However, 
presently, there are very few precise determinations of 
both the mass and radius of a compact object. There¬ 
fore, in the present paper we have proposed, and prelim¬ 
inary investigated, an alternative indirect method, that 
could help observationally distinguishing between differ¬ 
ent classes of compact objects, and their equations of 
state. This method is based on the information extracted 
from observations of the basic physical properties of mat¬ 
ter forming thin accretion disks around rapidly rotating 
neutron, quark and Bose-Einstein Condensate stars. 

Due to the presence of a strong gravitational field, 
all the astrophysical quantities related to the observable 
properties of the accretion disks, are dependent, and can 
be obtained from the metric of the central compact ob¬ 
ject 33]. Due to the major differences in the exterior 
space-time geometry, neutron, quark and Bose-Einstein 
condensate stars show, at least on the theoretical level, 
some very important distinct signatures with respect to 
the disk properties. Therefore, the observational proce¬ 
dure of the analysis of the electromagnetic radiation of 
accretion disks may allow to discriminate between neu¬ 
tron, quark and Bose-Einstein Condensate stars, by giv¬ 
ing some specific distinct signatures that could differenti¬ 


ate between compact objects described by different equa¬ 
tions of state. In the present paper we have obtained the 
physical parameters of the disk - effective potential, flux 
and emission spectrum profiles - for several equations of 
state of the neutron, quark and BEC matter, respectively. 

As one can see from the flux integral in Eq. (1M1) . as 
well as from the explicit expressions of the specific energy, 
specific angular momentum and angular velocity, given 
by Eqs. Ill]), (ESI) and (1271) . respectively, the rather dif¬ 
ferent characteristics of the radial flux distribution over 
the accretion disk, the disk spectra and the conversion 
efficiency are due to the important differences between 
the metric potentials of the neutron, quark and Bose- 
Einstein Condensate stars, respectively. Even if the total 
mass and the angular velocity are the same for each type 
of the rotating central object (neutron, quark or Bose- 
Einstein Condensate star), with the stars having similar 
values of D, E and L, the radiation properties of the ac¬ 
cretion disks around these compact general relativistic 
objects exhibit observable differences [33(. The physical 
reason for these differences is that the proper volume, 
and in turn the function y/—g, used in the calculation of 
the flux integral, is strongly dependent on the behavior of 
the metric component g rr = ( df/dr) 2 grr, and therefore 
on the geometry of the space-time. The latter expres¬ 
sion contains the derivatives with respect to the radial 
coordinate r of the metric components p(r) and 7 (r), 
respectively, via the coordinate transformation between 
the coordinates r and r, which are extremely sensitive 
to the slope of the functions p[r) and 7 ( 7 *) [33|. There¬ 
fore, although the inner edges of the disks are located 
at almost the same radii, the maximum amplitudes and 
the numerical values of the energy fluxes emerging form 
the disk surface, and propagating in any solid angle, may 
show considerable differences for different equations of 
state of the neutron, quark and Bose-Einstein Conden¬ 
sate matter. These essentially geometrical effects also 
give the distinctive features in the disk spectra for the 
various types of central stars. 

In our preliminary and idealized theoretical study of 
the accretion disk properties around Bose-Einstein Con¬ 
densate stars we have found a number of observational 
signatures distinguishing this class of stars from the neu¬ 
tron and quark stars. These specific properties are the 
distinct positions of the maxima of the flux, of the tem¬ 
perature distribution, and of the luminosity of the disk, 
the position of the inner edge of the disk, and the radi¬ 
ation efficiency conversion. In all three different classes 
of rotating stars we have analyzed these signatures do 
appear distinctly. Moreover, Bose-Einstein Condensate 
stars have a mass-radius relation that can also help in 
discriminating them with respect to other classes of neu¬ 
tron and quark stars. 

In the present paper, which represents a first step in 
the investigation of the complex astrophysical problem of 
the radiation emission from accretion disks around com¬ 
pact objects with different nuclear equations of state, we 
made the fundamental assumption that the inner edge 
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of the disk is located at the ISCO radius. This means 
that the ISCO position is only determined by the space- 
time geometry around the compact object. However, a 
number of important physical factors can affect the po¬ 
sition and location of the inner edge of the disk. For 
example, in the case of neutron stars, the inner edge of 
the disk is usually set by the magnetosphere around the 
star, and it is not determined by the metric of the space- 
time only. But, in order to simplify the discussion of 
the complex physics of the accretion disks, and by tak¬ 
ing into account the preliminary and purely theoretical 
nature of this study, in the present paper we consider 
only an ideal case, which allows us to start our analy¬ 
sis from very simple physical and astrophysical consid¬ 
erations. As already mentioned, an important physical 
parameter that strongly influences the ISCO position is 
the magnetic field of the compact central object, and its 
magnetosphere. The magnetic field lines originating from 
the central object can have a considerable effect on the 
evolution of the on the accretion processes, and on the 
disk structure. The effects of magnetic fields on the accre¬ 
tion processes in a compact, spherically symmetric geom¬ 
etry of Schwarzschild type (which describes the exterior 
geometry for both black holes and compact stars) were 
considered in [55J. The magnetic field was assumed to 
be asymptotically uniform, and axisynnnetric tidal struc¬ 
tures were also taken into account. Due to the presence 
of the magnetic field and of the tidal perturbations, the 
accretion disk shrinks in size, and the marginally sta¬ 
ble orbits shift towards the central object. On the other 
hand the presence of the magnetic perturbation leads to 
an increase in the disk radiation intensity from the accre¬ 
tion disk. However, the position of the maximum of the 
radiation does not change, and the radiation spectrum is 
blue-shifted. Interestingly, the conversion efficiency de¬ 
creases due to the presence of the magnetic fields and of 
the tidal perturbations. Alternative studies of the effects 
of the neutron star’s magnetic field on the exterior metric 
and the position of ISCO’s can be found in [56j . 

In the present paper we have proposed a method for 
discriminating between different types of compact objects 
that was proposed, and observationally tested, in the case 
of the black holes. Black holes have an event horizon (or a 
very special surface with similar properties) that prevent 
them to emit any form of electromagnetic radiation [57t | . 
This makes the radiation emission from the disk, and 
due to accretion processes, to be the dominant electro¬ 
magnetic energy emission mechanism, thus allowing the 
possibility of determining the central black hole proper¬ 
ties from the radiation spectrum of the disk, without any 
possible interference of the electromagnetic signals from 
the black hole. On the other hand, neutron stars have 
a solid radiation emitting surface (crust), which makes 
very difficult to distinguish between the radiation of the 
star and of the disk itself. That’s why the continuum- 
fitting method was not applied for the study of neutron 
stars [57}, and up to now there is no published obser¬ 
vational study of the possibility of testing the neutron 


star equation of state by using the electromagnetic spec¬ 
trum of a thin disk. On the other hand we would like to 
point out that the electromagnetic emissivity of the zero 
temperature gravitationally bounded Bose-Einstein Con¬ 
densate matter, in which all particles are in the quantum 
ground state, described by a single wave function, is very 
low [58]. This is because the effective plasma frequency 
of the condensed mater is very high. Therefore we may 
assume that the radiation emissivity of the BEC stars is 
very low, and they are very ’’black”. Consequently, the 
electromagnetic emissivity of the BEC star - disk system 
may be dominated by the disk emission. This situation is 
similar to the case of quark stars. Since quark matter has 
a very high plasma frequency w p , the photon emissivity 
of strange quark stars is very low |50[ . This is due to the 
fact that the propagation of electromagnetic waves hav¬ 
ing frequencies lower than lo p is exponentially damped. 
Hence, only photons produced just below a few fermi 
from the surface, with outwards pointing momenta can 
be emitted by the strange star. Hence the equilibrium 
photon emissivity from a strange star is negligible small, 
as compared to the black body one. Moreover, the spec¬ 
trum of the emitted equilibrium photons is very hard, 
with hu> > 20 MeV [59|. Of course BEC stars, as well 
as the quark stars, may have a crust (solid surface), rep¬ 
resenting a powerful source of electromagnetic radiation. 
The presence of such a crust would further complicate 
the possibility of discriminating between BEC stars and 
the other types of compact general relativistic objects. 

Another important point we would like to stress is that 
from an observational point of view not all the parame¬ 
ters of the BEC star-disk system can be determined inde¬ 
pendently from astrophysical observations. This means 
that a degeneracy between the equation of state and the 
nature of the compact object, and the values of some 
physical parameters of the model, could always exist. A 
full test of the BEC equation of state would require the 
precise knowledge of the mass, radius, and spin of the 
central object, as well as all the electromagnetic disk 
properties. In this ideal case, extremely difficult to be 
achieved from observational point of view, the fitting of 
the observational data could provide a convincing test 
of the nature of the central object. Probably in the near 
future no such increase in the precision of the astrophysi¬ 
cal observations will be achieved, thus making the direct 
determination of the EOS of the nuclear or condensed 
matter to be beyond the present observational capabil¬ 
ities. On the other hand one could expect an increase 
in the determination of the masses and radii for various 
compact objects, and these measurements may give some 
hints on the true nature of the equation of state of the 
dense matter. 

In conclusion, once the precision of the astrophysical 
data is drastically increased, the observational study of 
the thin accretion disks around rapidly rotating compact 
objects, and of their electromagnetic properties (flux, 
temperature distribution and luminosity), may provide a 
powerful tool in distinguishing between different classes 
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of dense stellar objects, as well as for discriminating be¬ 
tween the different equations of state of the dense matter. 
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